Irreversibility and Polymer Adsorption 
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FIG. 1: (a) Final irreversible layer structure. Chains high- 
lighted in bold: one belongs to the inner flattened layer (ujN 
surface contacts) the other to the outer layer (fN <C N con- 
tacts, loop size s ps n C ont//). (b) Late stage chain adsorption 
as surface approaches saturation and free supersites (clusters 
of n cont empty sites) become dilute. Chains cannot completely 
zip down. The minimum loop size s just connects two nearest 



neighbor supersites separated by l a 
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Physisorption or chemisorption from dilute polymer solutions often entails irreversible polymer- 
surface bonding. We present a theory of the non-equilibrium layers which result. While the density 
profile and loop distribution are the same as for equilibrium layers, the final layer comprises a tightly 
bound inner part plus an outer part whose chains make only fN surface contacts where N is chain 
length. The contact fractions / follow a broad distribution, P(f) ~ /~ 4 , in rather close agreement 
with strong physisorption experiments [H. M. Schneider et al, Langmuir 12, 994 (1996)]. 
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The validity of the laws of equilibrium statistical me- 
chanics hinges on ergodicity, the ability of a system to 
freely explore its phase space [Q. Many real processes, 
however, involve irreversible microscopic events such as 
strong physical or chemical bonding which invalidate 
ergodicity Equilibrium then becomes inaccessible and 
Boltzmann's entropy hypothesis is no longer applicable 
to calculate observables. Instead, the kinetics must be 
followed from their very beginning: the accessible region 
of phase space is progressively diminished as successive 
irreversible events freeze in an ever-increasing number 
of constraints. The state of the system at some time 
depends on the pocket of phase space to which it has 
become confined. 

The adsorption of high molecular weight polymers onto 
surfaces by its very nature frequently involves this kind 
of irreversibility (see fig. [I]). When an attractive surface 
contacts even a very dilute polymer solution there is a 
powerful tendency for dense polymer layers to develop 
|^|, [| because sticking energies per chain increase in pro- 
portion to the number of monomer units, N. This effect 
is exploited in many technologies such as coating, lubrica- 
tion, and adhesion. When the monomer sticking advan- 
tage e exceeds fcgT, available experimental evidence indi- 
cates that relaxation times become so large that the ph- 
ysisorption processes are effectively irreversible J| . This 
is a common situation. Many polymer species attach 
through strong hydrogen bonds || (e ~ 4fcsT) to silicon, 
glass or metal surfaces in their naturally oxidized states 
Pj , while DNA and proteins adhere tenaciously to a large 
variety of materials through hydrogen bonds, bare charge 
interactions or hydrophobic forces ||. In such situations 
layer structure is no longer determined by the laws of 
equilibrium statistical mechanics. The extreme exam- 
ple arises in chemisorption |?], |8| where covalent surface- 
polymer bonds develop irreversibly as in applications 
such as polymer-fiber welding in fiber-reinforced thermo- 
plastics and colloid stabilization by chemical grafting of 
polymers |^| . Generally, applications prefer the strongest 
and most enduring interfaces possible and irreversible ef- 
fects are probably the rule rather than the exception. 

Our aim in this letter is to understand the effect of 
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irreversibility on the structure of adsorbed polymer lay- 
ers (see fig. Q). Polymer adsorption phenomena are a 
major focus of polymer science, and though a few theo- 
retical and numerical works have addressed irreversibility 
H |l^, [n] the reversible case and the equilibrium layers 
which result are far better understood (5[ |J. Theory 
||, consistent with a number of experiments |jl2|| , pre- 
dicts each adsorbed chain in the equilibrium layer has se- 
quences of surface-bound monomers (trains) interspersed 
with portions extending away from the surface (tails and 
loops of size s). For good solvents the loop distribution 
Q(s) ~ s~ n / 5 and net layer density profile c(z) ~ z~ 4 ^ 3 
are universal. Equilibrium and ergodicity imply every 
chain is statistically identical. For example, for large N 
the fraction / of units which are surface-bound is the 
same for all chains to within small fluctuations and is 
no different to the overall bound fraction, / = Tbound/r. 
Here T is the total adsorbed polymer mass per unit area 
and Tbound the surface-bound part. 

How arc these universal features modified when the 
adsorption is irreversible? This question was explored 
in a series of ingenious experiments by the workers of 
refs. ^ who monitored polymethylmethacrylate (PMMA) 
adsorption from dilute solution onto oxidized silicon via 
hydrogen bonding with e w AksT. Measuring infrared 
absorption and dichroism, they monitored both T(t) and 
I\,ouiid (t) as they evolved in time and showed that early 




FIG. 2: (a) Predicted adsorbed polymer mass F versus 
surface-bound part Inbound- For chemisorption, Tbound ~ r 8//3 
initially, (b) Frequency histograms for fraction of bound mass, 
/. Experiment (grey) from ref. m Theory (empty), from pre- 
dicted distribution P(f) ~ /~ 4 ' 5 with / max < / ~ u, where 
values for / ma3C = 0.9 and u) = 0.47 were taken from ref. ^[ 



arriving chains had much higher / values than late ar- 
rivers and these / values were frozen in for ever. They 
modeled [|J this in terms of a picture where early arrivers 
lie flat and late arrivers having fewer available surface 
spots to adsorb onto are extended. The experimental 
/ values of the asymptotic layer followed a broad dis- 
tribution, shown in fig. ^](b). This succinctly quanti- 
fies the essential non-ergodic characteristic of these non- 
equilibrium layers: there are now infinitely many classes 
of chains, each class with its own particular statistics. 

In the following an initially empty surface contacting 
a dilute polymer solution with good solvent is consid- 
ered. We will calculate the kinetics of layer formation, 
Inbound (0 and r(t), and the distributions of / values and 
loop sizes in the evolving and final layer. The two cases 
of irreversible physisorption and chemisorption must be 
carefully distinguished. Define Q as the "reaction" rate 
between a monomer and the surface, given this monomer 
contacts the surface (see fig. ||). For physisorption, the 
attachment of a monomer is virtually instantaneous on 
reaching the surface so the effective value is diffusion- 
limited, Q w l/t a w 10 10 sec -1 typically, where t a is 
monomer relaxation time. Chemisorption processes are 



much slower, with typical values |]13| 10 
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sec . Consider a chain which, having diffused from bulk 
to surface, has just made its first attachment, i. e. just one 
monomer is irreversibly bonded to the surface (see fig. |^). 
We first treat the case of chemisorption, where the sub- 
sequent attachment of the remaining monomers is a pro- 
cess lasting seconds to hours and is thus experimentally 
accessible (all N monomers are assumed functionalized) . 

(1) Early stages: single chain adsorption and surface 
saturation. How does this chain adsorb down onto the 
surface? This depends on the exponent 9 governing the 
surface reaction rate k(s) for the s th monomer measured 
from the initial graft point (see fig. |^) 

k(s) « Q Z sur{ (s, N)/Z smi {N) « Q/s e , (s « N) . (1) 

Here Z sm f(N) and Z snT f(s, N) are the chain partition 
functions given one and two surface attachments, re- 
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FIG. 3: Chain adsorption commences with formation of an 
initial monomer-surface bond. For chemisorption, the reac- 
tion rate thereafter for the s th monomer from this graft point 
is k(s) ~ s~ e . Three modes of subsequent chain adsorption 
are theoretically possible: zipping (8 > 2); accelerated zip- 
ping, where occasional big loops nucleate new zipping centers 
(1 < 8 < 2); and uniform collapse (8 < 1). Chemisorption 
from dilute solution is accelerated zipping (8 — 8/5). 



spectively. Slow chemisorption allows sufficient time for 
chains to explore all configurations given the current con- 
straints frozen in by earlier reactions. Eq. (j|) states 
the reaction rate is proportional to the fraction of the 
grafted chain's configurations for which the s th monomer 
contacts the surface ||. Now in cases where 9 > 2, the 
total reaction rate 7?-totai ~ f-j* ds k(s) is dominated by s 
of order unity, i. e. monomers near to the first attached 
monomer will attach next. Thus, the chain zips down 
from the initial graft point. In contrast, for systems 
where 9 < 1 the upper limit dominates T^totai, i- e - a dis- 
tant monomer will react next; this implies a much more 
homogeneous chain collapse mechanism (see fig. |3|). 

The present situation is a self-avoiding polymer at a 
repulsive wall (we consider pure chemisorption, i. e. we 
assume a free energy advantage for solvent to contact 
the wall.) It turns out this case is intermediate between 
zipping and collapse. By relating 9 to other polymer ex- 
ponents at hard walls p3 we obtained the exact relation 



= 1 + 1/ 



(2) 



where v « 3/5 is the Flory exponent |l4| ] determining 
the polymer bulk coil size Rf — aN" in good solvent (a is 
monomer size). Thus 1 < 9 < 2 and 72-totai is dominated 
by its lower integration limit. We call this case accelerated 
zipping (see fig. ||). Zipping from the original graft point 
is accompanied by the occasional grafting of a distant 
monomer producing a loop of size s, say. This occurs 
after time r s w 1/ f s ds' k(s') « Q _1 s 3 / 5 . Each such 
new graft point nucleates further zipping, enhancing the 
effective zipping speed. Hence the entire chain adsorbs in 
a time ^adsorb = tjv ~ Q -1 ^ 3 / 5 , since by this time even 
the biggest loops have come down. Note this is much less 
than the pure zipping time « Q _1 N . Thus pure zipping 
must have been short-circuited by large loop adsorption 
events before it could have completed its course. 

During this accelerated zipping down, a characteris- 
tic (unnormalized) loop distribution f2t(s) develops and 
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the number of surface-bound monomers 7bound(i) grows 
from 1 to order N. We calculated these quantities by 
solving the detailed loop kinetics p|. These are rather 
complex, and here we present more accessible scaling ar- 
guments which reproduce the same results. Let us pos- 
tulate that after time t the only relevant loop scale is the 
largest to have come down, s max w (Qt) 5 / 3 , i.e. Q t (s) ~ 
(smax/s) a /smax for s < s max . Assuming a > 1, the to- 
tal number of loops L(t) « is dominated by small 
loops of order unity. Writing 7bound(*) = N^/tn) 13 we 
demand this be independent of TV for t <C tn (imagine 
sending the chain size to infinity; this would not affect 
the accelerated zipping propagating outwards from the 
initial graft point). This determines /3 — 5/3. Finally, 
since there are L(t) nucleating points for further zipping, 
ehWmdM ~ L, i.e. 7bound ~ Lt which fixes a = 7/5. 

We now sum over all chains which attached up to time 
t. The entropic disadvantage to touch the surface reduces 
the monomer volume fraction at the surface from the far 
field bulk value sur f = r<f> where the ratio of surface to 
bulk chain partition functions r = Z sur { (N) /.Zbuik (N) = 
1/N was calculated in ref. [lj. Then with 01 a 2 dT/dt = 
QN(j) surl and r bo und = 7bound(r/iV) we have 

T(t)a 2 = <j>Qt , r bound (i)a 2 = 4>N 3 / 5 {t/T N f/ 3 , (3) 

describing the early chemisorption layer for t < £ a dsorb = 
Q-ijV 3 / 5 . The loop structure of the partially collapsed 
chains is flt(s) ~ s~ 7 / 5 with maximum size s max = 
(Qi) 5 / 3 . This first phase may be long lived; e.g. for 
Q^ 1 = 1 sec, N — 10 3 then tn ps 20 mins. This becomes 
many hours for smaller Q values which are common. 

By time £ a dsorb zipping is complete and each chain is 
completely flattened onto the surface with fraction of ad- 
sorbed monomers f — lu. The species-dependent con- 
stant u is of order unity and reflects steric constraints 
preventing every monomer from actually touching the 
surface. In practice, we expect broadening of / val- 
ues about lu due to strong fluctuations, typical of mul- 
tiplicative random processes characterizing irreversibil- 
ity. For longer times each new chain zips down and 
rbound(i) = uT(t) with T given by eq. @. This pro- 
ceeds until ig at om wl/ {Q4>) when the surface is virtually 
saturated with a near-monolayer of flattened chains ||l7f. 

Consider now physisorption in its early stages. After 
attachment of its first monomer, the collapse of a single 
chain into a flattened structure now occurs as rapidly as 
monomers can diffuse a distance of order Hp, possibly 
accelerated by the attachments themselves. Thus we ex- 
pect the collapse time to be at least as small as the 
bulk coil relaxation time Tbuik (of order microseconds). 
Hence the collapse itself is probably experimentally un- 
observable, at least with the techniques of ref. |^. What is 
important is that in dilute solutions chains collapse into 
flattened configurations without hindrance from others. 
Moreover, we find that the probability a chain arriving 
from the bulk makes at least one bond before diffusing 



away is essentially unity even for a nearly-saturated sur- 
face. It follows that the attachment of chains is diffusion- 
controlled for essentially all times, a 2 r(i) w {(j)/ a){Dt) 1 / 2 
where D is center of gravity diffusivity. As for chemisorp- 
tion, Tbound = luT and adsorption produces a virtual 
monolayer of flattened chains. Surface saturation effects 
onset after time i s p a h t ys = T bulk (<£* / '4>f N 2 / 5 . 

(2) Late stages: the tenuously attached outer layer. 
Both chemisorption and physisorption processes fill the 
surface with completely collapsed chains, albeit in very 
different timescales t^ m and is at ys - By this stage the dis- 
tribution of surface-bound fractions is sharply peaked at 
f = lu. However, as saturation is approached free surface 
sites become scarce and late-arriving chains can no longer 
zip down completely. Suppose each chain-surface adhe- 
sion point consists in n con t attached monomers. The pre- 
cise value of n con t is sterically determined and is expected 
to be strongly species dependent. Then the surface den- 
sity of free "supersites" (unoccupied surface patches large 
enough to accommodate n cont monomers) is p S upcr ~ 
Ar bo und/"cont where Ar bo und = rg£ und - r bo und is the 
density of available surface sites and r b ^ iund is the asymp- 
totic density of bound monomers. Now as the sur- 
face approaches saturation so the density of supersites 
becomes small, p S upor <C l/(n. con ta 2 ), and their mean 
separation Z scp w p S upcr becomes so large that a late- 
arriving chain cannot find contiguous supersites to com- 
plete its accelerated zipping down. The minimum loop 
size s which can come down is that just large enough 
to connect two free supersites, i. e. as 3 ^ 5 — l scp whence 
s = (n cont /a 2 Arbound) 5 / 6 - Thus the final adsorbed state 
of chains arriving at this stage (see fig. [j](b)) consists of 
trains of ?i con t monomers separated by loops of order s 
units. For these chains <9 Arbound /c?Ar = / ~ n cont /s for 
large s, where Ar is the deviation from the asymptotic 
coverage T°° . Integrating this process up to saturation, 

a 2 Ar bo und = n cont (a 2 Ar/6) 6 , P(f) - Af^ (4) 

where / <C 1 and A is a constant of order unity [5j| . 
Adding this broad distribution of / values to the peak 
centered at / = to from the early stages gives the total 
distribution, shown in fig. |](b). It agrees rather closely 
with the experimental one of ref. |J shown in the same 
figure. The predicted rbound (r) profile (see fig. ||(a)) is 
also very close to the measured profile (i). 

Eq. (||) describes a tenuously attached outer layer 
(small / values) formed by late arriving chains, adding 
to the dense flattened layer formed at earlier times. The 
loop distribution of this diffuse outer layer is obtained 
from sVL{s)ds/Y ca = P(f)df whence 

fi(s) « a^s- 11 ' 5 , c{z)^z- i ' 3 (5) 

where the density profile followed from c = ilsds/dz 
evaluated at z = as 3 / 5 . 
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Finally, the kinetics of the total and bound coverages 
during the late stages are modified by saturation effects. 
For chemisorption the rate of attachment is directly pro- 
portional to the density of available surface sites, T ~ 

Ar bound ~ (Ar) 6 so Ar ~ t-vs an d Ar bound ~ t- & /\ 

In the physisorption case as discussed diffusion-control 
always pertains, V ~ i 1 / 2 , and thus the bound fraction 
saturates as AFbound ~ (1 — const. t/t^J? a ) & . 

In conclusion, we found that irreversible adsorption of 
polymer chains leads to final non-equilibrium layers ex- 
hibiting both similarities and profound differences com- 
pared to their equilibrium counterparts. The layer is a 
sum of a surface monolayer plus a diffuse outer part of 
thickness of order the bulk coil size with density profile 
c(z) ~ z~ 4 / 3 and loop size distribution O(s) ~ s~ n / D . 
Interestingly, these features are identical to those pre- 
dicted for equilibrium layers, including the precise expo- 
nent values. Prefactors are different, however, and we an- 
ticipate different values for physisorption and chemisorp- 
tion. To determine these necessitates accounting for 
topological constraints and fluctuations in empty surface 
site densities and other quantities, effects absent from our 
model. Note that although we did not explicitly treat ex- 
cluded volume interactions between an adsorbing chain 
and those previously adsorbed, we expect these to be 
unimportant because an empty site is correlated with a 
reduced surface loop density at that location. 

What is very different about irreversible layers is that 
individual chains in the layer are not statistically iden- 
tical: a given chain either belongs to the surface bound 
part and has order N surface contacts, or else the diffuse 
outer part. In the latter case the number of contacts, fN, 
is generally much less than N and its loop distribution is 
almost monodisperse with loop size s ~ 1//. In equilib- 
rium layers there is just one class of chain; parts of each 
chain lie bound to the surface, other parts extend into 
the outer layer and its loop distribution is the same as 
the layer's. In contrast, for irreversible layers there are an 
infinite number of classes, each with its own / value. The 
weighting for different values is universal, P(f) ~ /~ 4 / 5 
for small /. Practically, these differences have important 
implications for the physical properties of irreversible lay- 
ers; for example, the outer layer is much more fragile 
than the protected inner flattened layer. From a funda- 
mental point of view these systems provide a measurable 
example of how irreversible events progressively dimin- 
ish the available phase space volume and modify the en- 
tropy algorithm. For an equilibrium layer with e > kT 
this gives | F« F trans + F oam + F tra m for the free en- 
ergy. Here Ftrans = —ksT J dsil(s) ln[a 2 f2(s)] derives 
from loop translational entropy, F train = E tiain - TStrain 
is the contribution from trains and F osm is the osmotic 
part due to the solvent-swollen loops in the outer part of 
the brush. By comparison, for the non-equilibrium lay- 
ers both trains and loops are immobilized on the surface, 
ftrans = Strain = 0. The free energy is thus increased, 



F ~ F osm + E sur {. Its modified structure is expected, for 
example, to profoundly modify the interaction between 
polymer-covered surfaces as compared to the equilibrium 
case where the rearrangement of chains on the surfaces 
leads to characteristic force profiles ||. 
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